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Diode modeling for Schottky varactor frequency multipliers above 500 GHz is presented with special 
emphasis placed on simple models and fitted equations for rapid circuit design. Temperature- and doping-dependent 
mobility, resistivity and avalanche current multiplication and breakdown are presented. Next is a discussion of 
static junction current, including the effects of tunneling as well as thermionic emission. These results have been 
compared to detailed measurements made down to 80 K on diodes fabricated at JPL, followed by a discussion of the 
effect on multiplier efficiency. Finally, a simple model of current saturation in the undepleted active layer suitable 
for inclusion in harmonic balance simulators is derived. 

The research described in this publication was carried out at the Jet Propulsion Laboratory, California Institute of 
Technology, under a contract with the National Aeronautics and Space Administration. 


] Introduction. 

As diode multipliers are applied to higher frequencies (here used to mean above 500 GHz) the frequency 
dependent effects that can often be neglected at lower frequencies become more important. These include current 
saturation carrier inertia, and capacitance of the undepleted epitaxial region [1,2]. More accurate multiplier 
performance predictions benefit from including these effects in the models used. Furthermore, current and future 
space missions are being designed to operate the LO multiplier chains at lower temperatures, near 1 20 K for instance 
[3]. Again, for accurate modeling temperature effects must be considered. These include the effects on mobility, 
br^kdown' voltage, static MV current and current saturation. To aid in the rapid design of multipliers at many 
different power levels and frequency ranges, it is desirable to have models which are easy to implement using simple 
equations, with a minimum of arduous calculation, maximizing the speed of analysis. Thus, the more advanced 
models should still be fairly simple, suggesting the use of fitted equations. 
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Figure 1. Diode circuit model, (a) Simplified, (b) More complete 


As an illustration, Figure la 
shows a widely used diode model, 
which does a good job at lower 
frequencies, and can even work higher 
if well -calibrated. Figure lb indicates 
the more sophisticated model 
discussed in this paper. The elements 
of the model will be investigated in the 
following sections. In section H the 
temperature-dependent material 

properties are discussed: mobility, 

resistivity, avalanche current 
multiplication and breakdown voltage. 
Section III is an analysis of the static 
junction current including the effects 
of tunneling, and its effect on 
miiltinlipr pfTiripnrv In section IV a 


simple model for analysis of current saturation in the undepleted epi region is introduced. The conclusions are 
presented in Section V. 
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11 Diode Properties. 


A, Mobility. 

Schottky diode multiplier is signal absorption in the parasitic 
senes resist^ce depicted as R, m Figure la and and in Figure lb. TTie value oT these resXc“ is 
proporbonal to the resistivities in the active epitaxial layer, the highly-doped sub-layer, and ohmic contacts 
respwtively. Resistivity de^n^ on the low field mobility, which can either be calculated by Monte Carlo methods 
.""^ed p, Bl To facilitate calculation of the mobility for rapid design, a algebraic equation is 
dtsired, mcludmg the effects of temperature and donor concentration. The author is not aware of a good fit so the 
following expressions are suggested: ® ’ 
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[10]. Broken line: lOCK calculation, 
formula above are compared to the room temperature plot in [10]. 


where T is the temperature in Kelvins, no is the carrier 
concentration in cm ’, and 9 is the compensation ratio, defined 
[5] as the ratio of ionized acceptor to ionized donor 
concentrations. Since the normal donors in GaAs, such as Si, 
are shallow (~6 meV), as the carrier concentration increases the 
doncM- levels merge with the cmiduction band [9]. This occurs 
at a donor concentration of around 10** cm"’, and above that all 
donors are considered ionized except at extremely low 
tem{^tures, 10 K ot so. Since the range of temperatures 
considered here is above 50 K and dcmor concentrations are 
above 10 cm" , the earner concentration is assumed equal to 
effective donor concentration. No — N 4 , where No and are 
the donor and acceptor concentrations. 

The resistivity p is given by the usual formula: 

\/p = qpn^ ( 2 ) 

Figure 2 shows a plot of resistivity versus carrier concentration. 
Our room temperature and lOOK calculations based on the 


B. Avalanche current multiplication and breakdown. 

■ multiplication caused by impact ionization is impmtant in the reverse bias regime 

especially for the hi^ly doped diodes being used. This is the mechanism for junction breakdown, assuming the epi 
region IS thicker than the depletion width when it occurs. Avalanche multiplication increases the current injected 
one side of the depletion region by a factor M. Since impact ionization occurs at high fields, around 5x10’ to 
10 V/cm, It occurs first mdra- the anode at high reverse bias, that is, the current being multiplied is comprised of 
electrons penetrating the barrier from the anode contact. Under friese circumstances, M is given by [1 1]: 

1 ^ r 1 

^-T 7 = \f^n^^V-\(pC„-p^yx' dx ( 3 ) 

^0 L 0 

* 1 ^ field-dependent electron ionization coefficient, is the corresponding hole ionization coefficient, 
and IF is the depletion width at the given voltage. To determine the breakdown voltage, M is taken to be infinite, 


486 



12'^' International Symposium on Space Terahertz Technology^ 


and the right side of (3) is solved iteratively to determine the applied voltage and corresponding field distribution 
and depletion width that make it unity. 

Many researchers over the years have attempted to determine the ionization coefficients of GaAs [12-16]. 
For this work the measurements reported by Bulman et al. [14, 15] were used. BarafT [17] made calculations 
parameterizing the ionization coefficients at low temperature in terms of the optical phonon energy Er, the mean free 
path between phonon scattering events A, and the average ionization energy, These have been extended to all 
tempertures [18] using the assumption that E, varies with temperature as the band gap energy and that Er and A vary 
as; 


-lo U*''; 

where the zero subscripts refer to zero Kelvin values. The 0 K phonon energy is known to be 0.036 eV, so only 
values for £, and must be determined. The calculations in [17] have been fitted by Okuto and Crowell [18] to a 
simple formula. For our study the expression in [18] has been modified to fit BarafTs data more closely, which 
seems to model the measurements better. The modified equation is: 

aX- f ^ exp a -^la^ - x^ (5) 


Here a is either or pp as appropriate, /is a simple factor which accounts both for variations in measurments and 
for the unknown ratio of impact ionization cross-section to phonon emission cross-section [17]. The expressions for 
a, and jc are: 

/ ^ \1 266 / _ \ 


with E being the magnitude of the electric field. For this work values for the parameters that fit well to Bulman’s 
data are: = 2 eV and / = 0.44 for both electron and hole ionization coefficients, and Ao = 78.6 A forc^ and 71. 1 A 

for pp. As pointed out in [15, 19, 20] the ionization coefficients also depend on location. The effect of this was 
checked using the model of [15,20] , and it was found that for doping up to about 5xl0‘^ the difference in 
breakdown voltage with and without this correction is only a few tens of millivolts. 

Figure 3 shows the breakdown voltage 
calculated using the ionization coefficients from 
equation (5) in equation (6) at 300 K, compared to a 
plot from [21] based on PearsalFs coefficients [13]. 
Also shown as scattered points are measurements 
used to derive the coefficients in [14,15] as well as 
others as reported in [22, 23]. The fit is better to 
those using Bulman ’s data al lower doping, although 
the difference between the two calculations is very 
small for the more highly doped diodes. Also shown 
is a curve indicating the breakdown voltage at 100 K 
calculated as above. Again, the temperature 
dependence is small at the higher dopings. 
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Figure 3. Breakdown voltage cbping dependence. Scattered points 
represent measurements as described in text. 


through the Schottky barrier is the dominant static 
current mechanism in low temperature and in highly 
doped diodes. This is because at high doping levels 


the depletion width is narrow, giving the large number of electrons below the top of the barrier a substantial 
probablity of tunneling through. At lower donor concentrations the depletion width is great enough that the 
dominant current mechanism is thermionic emission over the barrier. Several approximate tunneling current 


calculations were developed by Padovani and Stratton and others [24-27], but due to the approximations in 
calculating the tunneling coefficient the results are not valid over all current regimes. Recently the transfer matrix 
technique has been used [28, 29] lo calculate the tunneling currents [30, 31]. This is a type of mode matching 
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algonthm wherein the barrier is divided up into a series of slices, each having a constant or linear potential profile 
After the transmissirai coefficient, 7i:£„) is calculated for each value of the nonnal componait of electron energy, 
the curroit density is calculated according to [32, equation (13)]. The transmission coefficient depends on the 
potential fa-ofile, including the barrier height, which is given for an anode at x = 0 by; 


\p = (t>b-qE^x\ 1 — — 


W j \67CE^£f^X 


with ErEo the semiconductor dielectric constant. The electric field under the anode is given ; 


whCTe V„ is the difference in potential between the Fami level and the bottom of the conduction band in the 
semiconductor bulk. The barrier height itself commonly appears to have a depoidoice on the applied voltage [10, 
33, 34]. For this work, the barrier height is characterized by an asymptotic “flat-band” height, 0fb, which is the 
barriCT height at zero electric field, when the conduction bands would be flat. Additionally, the b^er height is 
assumed to have a linear dependence on Emax‘’ 


0b =^FB ( 8 ) 

Since and are mutually dependent, they must be solved fw together as in [33], 



Figure 4. Block Diagram and photo of cryostat set up and wire-bonded diode. 


To check the validity and determine die effect of the tunneling current on multipliw paformance, sevia-al 
diodes fabricated at JPL have been measured over a wide temperature range in a specially designed cryostat [35]. 
The cryostat, shown in Figure 4, consists of a vacuum chamber «iclosing the device under test (DUT) cooled by a 
two-stage closed-cycle Helium refrigerator. The temperature of the DUT is adjustable to within ±1 K in the 35- 
325 K range using a temperature controller with two temperature sensors, a 25 W heata- and a bracket that 
mechanically connects the DUT to the 15 K second stage cold plate of the refiigeratOT. 

The HP 4155B Semiconductor Parameter Analyzer is connected to the DUT with four cables for a four- 
probe measurement, compensating for the series resistance of the lines (Figure 2). Shielded cables connect the 
analyzer to SMA adaptors at the vacuum flange of the cryostat where die signals are fed through into the chamber. 
Inside the cryostat short sections of flexible cable connect die feed-through to the ceramic test socket which is 
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thermally bonded through a copper block to the temperature-stabilized cold plate. The diode is epoxied to the test 
socket and connected with four 25 pm bond wires to the carrier as also shown in Figure 4. 

A diode fabricated a JPL on material doped to 3x10^^ with an anode size of 1.5x10.5 pm^ was measured 

from 100 to 300 K. The magnitude of 
the measured and calculated currents are 
compared in Figure 5. The current 
calculations included tunneling and 
avalanche multiplication as described 
earlier. Values of 0fb = 0-92 eV and a 
= 30 A in equation ( 11 ) were used to fit 
the measurements. The fit is pretty 
good, but it is desirable to have 
independent justification for these 
values. The flat-band barrier height is 
dependent on the metal/semiconductor 
system in the contact. Using forward 
bias measurements of the reverse 
saturation current and ideality factor, the 
flat-band barrier height and effective 
barrier height, ^ ^ determined as 

discussed in [36-39] based on the fact 
that as the applied voltage approaches the built-in voltage (flat-band), all current mechanisms except thermionic 
approach zero. This includes tunneling (because the bulk of electrons are near the barrier top in energy) and any 
models described by a linear field dependence (since the field approaches zero). This type of analysis was also 
performed on the above diode, with the results shown in Figure 6 a. Since the JPL anodes are fabricated from 
Ti/Pt/Au, the flat-band barrier height should be compared to those metals. In [38] the reported heights are around 
0.83 eV for Ti, 0.92 eV for Au and 0.99 eV for Pt. The measured values fit a line 0fB = * • •-•37' meV, T being the 
temperature in Kelvin. It is possible that some intermixing may be taking place between the Ti and the Pt, but this 
cannot be said for certain. As to the a== 30 A factor, many theories have been proposed explaining non-zero a as 
due to the imperfect nature of the metal-semiconductor interface [33, 34, 40-44]. Again, there is not enough data to 
pin this down. 



Temperature (K) Temperature (K) 


Figure 6 . Characteristics of 3el7 diode, (a) Barrier height and ideality. The solid line is a model 
C>p 0 = 0.95-0.0003T eV.( b) Series and shunt resistances. 

The series resistance and shunt resistance can also be found from the I/V data. These are plotted in Figure 
6 b. As expected, the series resistance R, is independent of temperature (except for deviations ascribable to 
measurement uncertainty). The shunt resistance, Rp, however rises sharply at low temperatures. This is consistent 
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with a model wherein the shunt conduction is a surface “hqjping” phenomenon, where the carriers do freeze out as 
the tempCTature drops. Generation-recombination currait was considered, but it is propOTtional to the product of 
depletion width IF and intrinsic carrier concentraticm n,. For GaAs in general «, is very low and for highly doped 
diodes W is small, so this current mechanism is insignificant in this case. 


Reverse; F«Tward: 



Voltage (V) 

Figure 7. Fitted diode models and measuiements, 

100 4& 300 K, Markers indicate fit models. 



Figure 8. Efficiency of doubler with and without 
including tunneling current. 


To evaluate the effect of these departures from thermionic behavior on multiplier efficiency, a simple 
model consisting of two diodes back-to-back was incOTporated into a commercial harmonic balance simulator. A 
schematic is shown in Figure 7. One diode models the forward current, the otha- the reverse. The fit of the modeled 
current magnitude to the measurements is also shown in Figure 7. The fit is good for the forward bias, and not so 
good for reverse. However, the fit at the higher reverse currents is good, and the lower curraits should have little 
effect on the performance of the multiplier. 

Performance calculations of a 200 GHz multiplier with and without the reverse tunneling component are 
shown in Figure 8. As can be seen, at room temperature there is a significant reduction in efficiency when tunneling 
current is considered, but at 100 K there is no substantial difference. There is also no observed difference when the 
same calculation is performed at 800 GHz. These can be explained by noting that all Schottky multipliers suffer 
from two loss mechanisms: series resistance loss, vAiich dominates, and conduction current loss. As the temperature 
IS lowered, the conduction currait drops (at a ^ven voltage) so the conduction current contribution to the loss is 
greatly reduced. Hence, little difference in efficiency occurs with or without tunneling at low temperatures. At 800 
GHz the efficiency is low and also greatly dominated by the series resistance loss, so the exact conduction currait 
mechanism has little effect on efficiency. 


IV Current Saturation. 


As the electric field in a GaAs sample is increased, the velocity and hence the current reaches a peak, tiien 
begins to decline due to electrois fining enough enagy from the field to scatta into the upper, low-mobility 
valleys. This acts as a current limit in the undepleted epitaxial region of the diode. Current saturation reduces the 
efficiency of the diode as a multiplier because it decreases the ability of the charge at the edge of the depletion 
region to move around at the signal frequency and modulate the capacitance, goiaating the varactor non-linearity. 
There are also time constants associated with the transfer, and at the high frequencies being considaed hae they 
should be taken into account. 

There are several equations used to describe the static current saturation profile. The simplest and oldest is 
the Kramer and Mircea formula [47]: 


l + (F/F„X 


( 12 ) 


where u is the velocity, F the magnitude of the electric field, the ultimate saturation velocity at about 20 kV/cm, 
and Fn is a charactaistic field which determines whae the peak velocity occurs. Usually current saturation,’ 
including transient effects have been modeled using Monte Carlo simulations (see, for example [45, 46]). Since 
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Monte Carlo calculations take large amounts of computer time, it is desirable to incorporate these effects into a 
harmonic balance (HB) circuit simulator using a simple model which can be integrated using die Runge-Kutta t;^e 
integrators normally used in HB simulators. We propose to use a time constant based formulation somewhat similar 
to that introduced in [48]. The epi current is divided between two resistances representing the dominant two 
conduction band valleys in which electrons travel, as illustrated in Figure 9. Defining po as the lower valley 
mobility, (about 4000 cm^/Vs or higher at reduced temperatures), and /// is the upper valley mobility (about 400 
cm^A's) and no and n, as the corresponding valley populations, the velocity in equation (12) can be written: 


u{f)= 


( 




V 




(13) 



Figure 9. Simplified model of current saturation. 


with n the total electron concentration in the undepleted epi. 
Then, several coupled differential equations are used to 
represent the time-dependent behavior of the velocity. The 
upper valley population is described by: 

^ «i.(^)~»i(0 ( 14 ) 

dt T 

where ni,(V) represents the static population of the upper valley, 
derived by combining equations (12) and (13) and noting that 
the total electron concentration is the sum of the populations of 
the two valleys, i.e. « = + rtj. The time constant, r is 

dependent on ihe voltage across the epi, Vtot^ 

The current through the inductor representing carrier 
inertia has the usual equation: 

di „ 5 , 

dt L 


where L, is the carrier inertia inductance, equal to: 


' q^nA 

with m* the effective mass, w the thickness of the undepleted epi, and A the area of the diode. The total current 
through the inductance is simply: 

i= nuA. 07 ) 

To simulate the full diode, of course, other equations must be included to model the epi capacitance, as well as the 
junction capacitance and conducticm current. 

As a check on the validity of the model, a test on a sample “undepleted epi” was made. It was assumed that 
this sample was of constant thickness, so that the equations above could be used, uncoupled from the epi 
capacitance, since the its current would depend on V,oton\y, which is given as a bound^ condition. In that case, 
equations (13) through (17) can be integrated using a numerical integrator with any given waveform for Vtot- In 
Figure 10 these calculations are compared with those in [46] for steps from zero electric field to the indicated values. 



Figure 10. New current saturation model transient 
response (smooth curves) and results 
from [6] (rough curves). 


Figure 1 1 . Sinusoidally pumped current saturation 
model (smooth line) compared to results 
from [46] (rough) and resistive approx- 
imation (sine wave). 
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It^ found *at values for r ranging linearly from 0.2 ps at the 20 kV/cm to 3 ps at 1 kV/cm gave a good fit. As a 
further test, the^e model is compared in Figure 1 1 with a sinusoidal 60 GHz excitation of 10 kV/cm amplitude. 
Again, Ae match is good. Also plotted is the current through a resistance having a mobility one-third that of GaAs 
(hence three times the resistance), which gives an estimate of the magnitude in performance reduction 

't appears that the above model, simple as it is, is impossible to incorporate into commercial HB 
simulators. Work is continuing in this direction, and a version of Ae Siegel/Kerr reflection algoriAm HB simulator 
I4vj IS in progress so this model can be used to actually design multipliers. 

V Conclusions. 

Several suggestions have been made for improving Ae modelmg of Schottky diodes for high doping and 
frequency, mid low temperatures. The effects exammed have been low-field mobility, breakdown, conduction 
cuirent, ^er hei^t and current saturation. Most of Ae effects serve to merely reduce Ae achievable efficiency 
but not chMge Ae designs Aemselves very much, specifically Ae matchmg impedances. This is because the mosi 
important design parameter is Ae reactance tunmg, dominated by Ae Aode capacitance, which has been accurately 
model^ previously. However, the current saturation serves to greatly increase Ae effective series resistance which 
IS Ae domin^t loss mechanism. It is possible to estimate Ae magniAde of Ae effect by using a series resistance 
much larger Aan Ae measured DC values m Ae performance estimates [50]. However, by calculatmg Ae behavior 
of current sat^tion wiA variations in dopmg and temperature, it should soon be possible to improve the 
optimization of the doping and anode sizes of high frequency multiplier diodes. 
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